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Abstract 

We report a Chandra high-energy grating detection of a narrow, redshifted absorption line superimposed on 
the red wing of a broad Fe K line in the z — 0.297 quasar E 1821+643. The absorption line is detected at a 
confidence level, estimated by two different methods, in the range ~ 2 — 3cr. Although the detection significance 
is not high enough to exclude a non-astrophysical origin, accounting for the absorption feature when modeling 
the X-ray spectrum implies that the Fe-K emission line is broad, and consistent with an origin in a relativistic 
accretion disk. Ignoring the apparent absorption feature leads to the conclusion that the Fe-K emission line is 
narrower, and also affects the inferred peak energy of the line (and hence the inferred ionization state of Fe). If 
the absorption line (at ~ 6.2 keV in the quasar frame) is real, we argue that it could be due to gravitationally 
redshifted Fe xxv or Fe xxvi resonance absorption within ~ 10 — 20 gravitational radii of the putative central 
black hole. The absorption line is not detected in earlier ASCA and Chandra low-energy grating observations, but 
the absorption line is not unequivocally ruled out by these data. The Chandra high-energy grating Fe-K emission 
line is consistent with an origin predominantly in Fe i-xvn or so. In an ASCA observation eight years earlier, the 
Fe-K line peaked at ~ 6.6 keV, closer to the energies of He-like Fe triplet lines. Further, in a Chandra low-energy 
grating observation the Fe-K line profile was double-peaked, one peak corresponding to Fe i-xvn or so, the other 
peak to Fe xxvi Lya. Such a wide range in ionization state of Fe is not ruled out by the HEG and ASCA data 
either, and is suggestive of a complex structure for the line-emitter. 
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1. INTRODUCTION 

The high-luminosity (L2_io keV ~ 3 — 4 x 
10 45 ergs s" 1 ) quasar E 1821+643 (z = 0.297) 
is one of the highest redshift quasars to exhibit 
strong Fe Ka line emission. The source has been 
observed by every X-ray astronomy mission since 
EXOSAT, but until recently, the contribution to 
the Fe Ka line emission from the cluster that the 
quasar is located in, was highly uncertain (e.g. 
see Saxton et al. 1997 and references therein). 
Fang et al. (2002) finally showed, from a Chan- 
dra observation, that the cluster makes a negligi- 
ble contribution to line emission at 6.4 keV (the 
energy of Fe I Ka) , and at most ~ 3% at Fe xxvi 
Lya (~ 6.9 keV). 

In this paper we show that one interpretation 
of the Chandra data for E 1821+643 is that the 
Fe Ka emission line is relativistically broadened, 
with an absorption feature at ~ 6.2 keV (quasar 
frame) superimposed on the broad emission line 
profile. Whilst blueshifted absorption features 



are common in active galactic nuclei (AGNs), a 
redshifted absorption feature has been reported 
for only one case so far (NGC 3516, Nandra et 
al. 1999). The absorption feature in E 1821+643 
could be due to resonance absorption in highly 
ionized Fe, in the form of an inflow, matter cross- 
ing the line-of-sight obliquely, or an outflow that 
is close enough to the putative central black hole 
that the absorption line suffers a strong gravi- 
tational redshift. The latter scenario would of 
course have very important implications for the 
study of the central engine with future missions. 
Whatever the origin of the absorption line, its 
presence in the HEG spectrum affects modeling 
of the Fe-K emission line at ~ 6.4 keV and there- 
fore cannot be ignored. The absorption may be 
of a transitory nature since it is not detected in 
another Chandra observation, nor was it detected 
with ASCA. However, each of these latter data 
sets were not as sensitive to the absorption line 
and we will show that they do not unequivocally 
rule it out. Our analysis also shows a significant 
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detection of Fe-K line emission at ~ 6.9 keV in 
Chandra low-energy grating (LEG) data (in ad- 
dition to the line at ~ 6.4 keV). Likely to be 
associated with Fe xxvi Lya, this line is weak 
in the HEG data and not detected in the ASCA 
data, but both data sets are statistically consis- 
tent with the LEG data. The overall line profile 
in the ASCA data has a broad peak at ~ 6.6 keV 
but multiple emission lines cannot be resolved. 

The paper is organized as follows. In §2 we de- 
scribe the Chandra high-energy grating data and 
its reduction; in §3 we describe modeling of the 
Fe-K emission and absorption lines in these data, 
including a detailed discussion of the statistical 
significance of the absorption line. In §4 we dis- 
cuss the Fe-K emission and absorption features 
with respect to two other important data sets for 
E 1821+643, namely from observations with the 
Chandra LEG and with ASCA. In §5 we discuss 
possibilities for the origin of the absorption line, 
and finally, in §6 we present our conclusions. 

2. THE Chandra HEG DATA 

E 1821+643 was observed with Chandra start- 
ing 9 February, 2001, from UT 14:28:55, for a du- 
ration of ~ 100 ks. This Chandra observation was 
made with the High-Energy Transmission Grat- 
ing (or HETGS - Markert, et al. 1995) in the 
focal plane of the High Resolution Mirror Assem- 
bly. The Chandra HETGS affords the best spec- 
tral resolution in the ^6 — 7 keV Fe-K band cur- 
rently available (~ 39 eV, or 1860 km s" 1 FWHM 
at 6.4 keV). HETGS consists of two grating as- 
semblies, a High-Energy Grating (HEG) and a 
Medium-Energy Grating (MEG), and it is the 
HEG that achieves this spectral resolution. The 
MEG spectral resolution is only half as good as 
that of the HEG. The HEG also has a higher ef- 
fective area in the Fe-K band, so our study will 
focus principally on the HEG data. The HEG 
and MEG energy bands are 0.4-10 keV and 0.7- 
10 keV respectively, but the effective area falls 
off rapidly with energy near both ends of each 
bandpass. The Chandra data were processed with 
version R4CU5UPD14.1 of the processing soft- 
ware, CALDB version 2.2 was used, and the tele- 
scope responses made using ciao 2.1.3 1 . Oth- 
erwise, HEG and MEG lightcurves and spectra 
were made exactly as described in Yaqoob et al. 
(2003a). We used only the first orders of the 

1 http://asc.harvard.edu/ciao 



grating data (combining the positive and nega- 
tive arms, unless otherwise stated). The zeroth 
order data were piled up and the higher orders 
contain much fewer counts than the first order. 
The mean count rates (in the full energy band of 
each grating), over the entire Chandra observa- 
tion were 0.1346 + 0.0023 ct/s and 0.2938 + 0.0027 
ct/s for HEG and MEG respectively. HEG and 
MEG spectra extracted over the entire observa- 
tion, resulted in a net exposure time of 99.620 
ks. Background was not subtracted since it is 
negligible in the energy region of interest. The 
source flux showed little variability over the en- 
tire duration of the campaign. For example, for 
HEG plus MEG lightcurves binned at 1024 s, the 
excess variance above the expectation for Pois- 
son noise (e.g. sec Turner et al. 1999) was 
(—0.4 ± 3.5) x 10~ 4 , consistent with zero. 

3. SPECTRAL FITTING OF THE 
BROAD Fe-K EMISSION LINE AND 
NARROW ABSORPTION LINE 

We used XSPEC vll.2 for spectral fitting (Ar- 
naud 1996), omitting the 2-2.5 keV band due 
to systematic residuals around sharp changes 
in the effective area of the telescope - see 
Yaqoob et al. (2003a). The response ma- 
trices acisheglD1999-07-22rmfN0004.f its and 
acismeglD1999-07-22rmfN0004.f its (for the 
HEG and MEG respectively), combined with the 
telescope response files described above, were 
used to fold models through the instrument re- 
sponse and thereby directly compare predicted 
and observed counts spectra. Therefore, quoted 
fit parameters do not need to be corrected for 
instrumental response. For fitting purposes, the 
HEG and MEG spectra were binned at O.OlA 
and 0.02A respectively. The bin size is compa- 
rable to the FWHM spectral resolution (0.012A 
and 0.023A for HEG and MEG respectively). For 
clarity, spectral plots in this paper show larger 
bin sizes than were used in the fitting (the bin 
size will be given case by case). 

Since some energy bins may contain zero or few 
counts, the C-statistic was used for minimization 
and, unless otherwise stated, all statistical errors 
correspond to 90% confidence for one interesting 
parameter (AC = 2.706). By definition, calcula- 
tion of the C-statistic requires only knowledge of 
the number of counts in a bin, but for spectral 
plots, the error bars shown correspond to asym- 
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Figure 1 . The HEG spectrum of E 1821+643 (corrected for instrumental response), binned at 0.04A. The source redshift is 
0.297. The solid line is the best-fitting power-law model fitted over the observed energy range 1.0-2.0, 2.5-9 keV (omitting 
the 2.0-2.5 keV region where the mirror response is not well modeled). The energy range 1-9 keV corresponds to 1.3- 
11.7 keV in the rest-frame. The inset shows the ratio of the same spectrum to the above power-law model in the Fe K 
region. An asymmetric, broad Fe-K line is apparent, with a deep absorption feature superimposed on the profile. The 
vertical dotted lines correspond (from left to right) to the energies of the following transitions: Fe I Ka, Fe XXV Is 2 — ls2p 
forbidden and resonance lines, and Fe XXVI Lyct (which correspond to rest-frame energies of 6.400, 6.636, 6.700, and 
6.966 keV respectively). 



metric errors calculated using the approximations 
of Gehrels (1986). All model parameters will be 
referred to the source frame using z — 0.297, un- 
less otherwise stated. 

We fitted the HEG spectrum in the 1-9 keV 
band, in the observed frame (except for the 2- 
2.5 keV band). Galactic absorption was not in- 
cluded in the model fitting since it has a negligible 
effect above 1 keV. First we fitted a simple power- 
law model. This model is shown in Fig. 1 over- 
laid on HEG data that have been corrected for 
the instrument resolution and binned at 0.04 A. 
The inset in Fig. 1 shows the ratio of data to 
model. It can be seen that there is broad, asym- 
metric excess emission centered around ~ 4.9 keV 
(<~ 6.4 keV quasar frame), likely corresponding 



to a low-ionization state Fe Ka line. A similar, 
but narrower emission line has been reported by 
Fang et al. (2002) using the same data. How- 
ever, there also appears to be a deep absorption 
feature, up to about 0.04A wide (~ 75 eV at the 
observed energy), superimposed on the red wing 
of the emission line, at ~ 4.8 keV in the observed 
frame (<~ 6.2 keV in the quasar frame). The ex- 
act value of the centroid energy of the feature, 
and associated statistical errors, depend on the 
model of the emission line and these matters are 
discussed in §3.2 and §3.3. 

3.1. Reality of the Absorption Line 

We investigated whether the apparent absorp- 
tion line could be due to an instrumental arti- 
fact or a statistical fluctuation. Wc will describe 
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Figure 2. Detection of the high-energy absorption line in 
the separate +1 (green crosses) and —1 (red open circles) 
orders of the Chandra HEG, compared with the spectrum 
summed over both grating arms (black filled circles). The 
absorption line is probably due to redshifted He-like or 
H-like Fe absorption, equivalent to an inflow velocity of 
~ 2 X 10 4 km s -1 or ~ 3 X 10 4 km s — 1 respectively (see 
text). Shown are the ratios of spectral data (binned at 
0.04A) to the best-fitting power law fitted to the spectrum 
combined from both +1 and —1 orders of the grating (i.e. 
the same power-law model as shown in Fig. 1). 



this analysis in considerable detail since there 
are important physical consequences if the ab- 
sorption line is real. Shown in Fig. 2 are the 
data/model ratios for the separate +1 and —1 
orders, compared with the combined first-order 
data/model ratio (at bin sizes of 0.04A), where 
the model is the simple power law fitted to the 
combined first-order data, as shown in Fig. 1 and 
described above. Fig. 2 shows that the absorp- 
tion line is present, at some level, in both the 
negative and positive arms of the grating. Since 
the physical distance on the CCD detector of a 
photon event from the zeroth-order image deter- 
mines the energy of that photon event, and the 
energy bin containing the absorption line in the 
negative and positive arms have completely dif- 
ferent physical locations on the CCD array, an in- 
strumental artifact cannot explain the absorption 
line. Also, the effective area of the HEG is smooth 
over the region the absorption line (centered on 
~ 4.8 keV, observed frame) is detected: the com- 
bined +1 and —1 order effective area varies by 



no more than 6% over the observed energy range 
4.6-5.0 keV. Yet the photon number in the 0.04A 
energy bin in which the line is observed, deviates 
by at least 40% compared to adjacent bins. A 
mis-calibration of the effective area would have to 
be present on extremely localized regions of dif- 
ferent CCDs, corresponding to precisely the same 
energy bin. There arc no other instrumental ef- 
fects that we are aware of that could explain the 
absorption feature cither. 

To quantify the possibility that a statistical 
fluctuation occurred in both arms of the HEG, 
in the same energy bin, to conspire to yield an 
apparent absorption line like that observed, we 
employed two different methods to assess the sta- 
tistical significance. The question we are seeking 
an answer to is the same in both cases. Namely, 
what is the probability of obtaining a statisti- 
cal deviation in photon counts as large as that 
observed for the absorption line in one arm of 
the grating in any energy bin in the energy range 
searched, and then to observe a deviation as large 
as that measured in the other arm of the grating, 
this time in the same energy bin? 

The first method we used is 
model-independent. For each first-order arm of 
the HEG, we calculated the fractional deviation 
of the counts in the 0.04A bin containing the ab- 
sorption line, relative to the mean counts in the 
two bins either side of it. We must account for 
the total number of bins that were searched that 
resulted in the detection of the absorption fea- 
ture, given that its energy does not correspond 
to a known atomic transition. For this, we used 
a very conservative assumption, counting all the 
energy bins that were used in spectral fitting. So, 
in the 1-2 keV and 2.5-9 keV range, for the +1 
order, we calculated, for every bin, the Poisson 
probability of obtaining the same fractional devi- 
ation, say X + (= 58.1%), as for the absorption 
line in that arm. For this, we assumed the mean 
expected counts to be equal to the average of the 
two bins surrounding the target bin. Then, by 
combining these probabilities, we calculated the 
total probability of observing the deviation, X + , 
in any of the energy bins in the energy range 1— 
2, 2.5-9 keV. Next, we calculated the probability 
in the other HEG arm (—1 order) of obtaining 
the observed deviation, A„ (= 47.9%), in the ac- 
tual absorption-line bin, relative to the average 
counts of the two surrounding bins. We obtained 
a final probability of obtaining the actual devia- 
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tions that were observed, X + and X_, in the +1 
and —1 arms of the HEG respectively (and simul- 
taneously), of 2.38 x 1CT 3 (i.e. this is the proba- 
bility that both detections were due to statistical 
fluctuations). This corresponds to an overall sig- 
nificance of 3.04cr for detection of the absorption 
line. 

We went further and estimated a significance 
using an even more conservative assumption. 
This was based on the possibility that both bins 
surrounding the absorption-line bin, in each HEG 
arm, could both have been higher than the expec- 
tation values, due to statistical fluctuations, thus 
artificially enhancing the depth of the absorption 
line. Therefore, we repeated the above calcula- 
tions after reducing the counts in all bins by la, in 
both HEG arms, except for the counts in the ab- 
sorption line bins. Here, we used 1 + ^/(0.75 + N) 
for the lcr error on the counts, N (see Gehrels 
1986). Strictly speaking, this is not correct for 
small N, but the bins either side of the absorp- 
tion line in each arm of the HEG have between 34 
and 43 counts, so the effect of the approximation 
on the final result is unimportant. (Note that the 
MEG has between 17 and 29 counts per bin in 
the bins either side of the absorption line, in each 
arm). Using this more pessimistic prescription 
we get 1.95cr for the significance of the absorp- 
tion line, which is consistent with the fact that 
we have effectively artificially reduced the level of 
the reference 'continuum' by ~ lcr. 

The second method we used to assess the sig- 
nificance of the line employed Monte Carlo simu- 
lations. We took a model for the continuum and 
emission line, and folded it through the instru- 
ment response, so that we could measure devia- 
tions of the data relative to this model. We used 
the power-law plus Gaussian emission line model 
described in §3.2 and Table 1, and set the absorp- 
tion line equivalent width (EW) in that model to 
zero. Obviously the probabilities calculated from 
the Monte Carlo simulations will depend some- 
what on the model but we will see that the re- 
sults are consistent with the model-independent 
method described above. 

We measured negative deviations relative to 
the continuum plus emission-line model of 20.5% 
and 38.2% in photon counts, in the 0.04A energy 
bin containing the absorption line, in the +1 and 
— 1 HEG orders respectively and used these as in- 
puts to the Monte Carlo calculations. We ran 10 6 
simulations for every energy bin, for each HEG 



Table 1 

Gaussian Line Fits to Chandra HEG Data for 
E 1821+643 



Parameter 


Measurement 


C-statistic 


1017.0 


Degrees of freedom 


968 


Emission Line Center Energy (keV) 


n oi +0.06 

U-Ol_ Q 5 


Emission Line Width, a (keV) 


0.18±S'o4 


Emission Line Velocity 


20145±Sgg 


(FWHM, km s" 1 ) 


5.2±i; 3 2 


Emission Line Intensity 


(10~ 5 photons cm~ 2 s _1 ) 


296t£ 


Emission Line EW (eV) 


Absorption Line Center Energy 


(. OOQ + 0.011 


(keV) 




Absorption Line Gaussian Width 


027+ ' 011 

U - UZ ' -0.007 


7 ( keV) 


3040l 7 8o° 


Absorption Line Velocity Width 


(FWHM, has- 1 ) 




Absorption Line Equivalent Width, 


EW (cV) 




Power-Law Photon Index T 


i po+0.03 
1 -° z -0.03 



Simple power-law model plus a Gaussian emission 
line, and a Gaussian absorption line, fitted to the 
Chandra HEG data (see §3.2 for details). All pa- 
rameters (except 2-10 keV flux) are referred to the 
quasar frame (z = 0.297). Errors are 90% confidence 
for one parameter (AG = 2.706). Velocities have 
been rounded to the nearest 5 km s _1 . 



first-order arm. We used the most conservative 
assumption again, for the number of energy bins 
searched, namely the 1-2 and 2.5-9 keV range. 
We measured the deviations relative to the model, 
in each bin, and in each simulation, and thus ob- 
tained a final probability that the observed ab- 
sorption line is due to chance statistical fluctua- 
tions, of 2.44 x 10~ 2 , or a significance of 2.25(7. 
This is consistent with the estimate of 2-3 a from 
the model-independent method described above. 
We note that these are all very conservative esti- 
mates of the significance, because in principle, we 
could have found the absorption line in the real 
data by examining a smaller number of energy 
bins. 
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From Fig. 1, it might appear as if there are de- 
viations as strong, or stronger, than those for the 
reported absorption line, notably at ~ 2.5 keV, 
and above ~ 6 keV (both observed frame). We 
remind the reader that 2-2.5 keV is the region 
of greatest changes in effective area, mainly due 
to the X-ray telescope (XRT), so even small in- 
accuracies in the calibration of the effective area 
can result in apparent emission or absorption fea- 
tures. Indeed, one finds that the deviations at 
~ 2.5 keV are not consistent in the +1 and —1 
HEG arms. In fact, we already mentioned in §3 
that we did not use the 2-2.5 keV data in any 
of our analysis or simulations, due to the difficult 
calibration of this region of the XRT response. 

On the other hand, the deviations above ~ 
6 keV in the +1 and —1 HEG orders do corre- 
late with each other and may represent real struc- 
ture. For example, Fe-K edges are in the right 
energy range (6 keV corresponds to ~ 7 keV in 
the quasar frame). However, the effective area of 
the HEG (and MEG) drops quickly above 6 keV, 
and the statistical significance of these deviations 
is much less than that of the reported absorption 
line, and is not sufficient to warrant further in- 
vestigation. 

We also investigated whether the absorption 
line is detected in the MEG data. A direct com- 
parison of the MEG and HEG data is given in 
Fig. 3. Shown are the ratios of MEG and HEG 
data to a simple power-law continuum fitted in 
the 1.0-2.0, 2.5-9 keV bands. Although no ab- 
sorption line is evident in the MEG data, it is 
worth remembering that the MEG effective area 
is a factor 1.5 less than that of the HEG at the 
energy at which the HEG absorption line is de- 
tected. We quantified the probability that the 
MEG would not detect the absorption line, given 
the HEG measurements. The observed number 
of counts in the 0.04A bin that would contain the 
absorption line in the MEG spectrum is 43. Tak- 
ing the best-fitting HEG power-law plus Gaus- 
sian emission line, and absorption-line model (see 
§3.2) and folding it through the MEG instrumen- 
tal response function, and fitting for the overall 
cross-normalization difference between HEG and 
MEG (<~ 8%), we estimated the mean expected 
counts in the absorption-line bin to be 32.2. The 
difference between expected and measured counts 
in that bin is only ~ 1.1a. The Poisson proba- 
bility of actually obtaining 33 or more counts in 
that bin is then 46.6%. In other words, this is the 
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Figure 3. A comparison of the HEG data (black) with 
the MEG data (red) in the Fe-K region in E 1821+643. 
Shown are the ratios of spectral data (binned at 0.04 A) 
to the best-fitting power-law fitted to the HEG spectrum 
combined from both +1 and —1 orders of the HEG (i.e. 
the same power-law model as shown in Fig. 1 and Fig. 2). 
Although the absorption line does not appear to be sig- 
nificantly detected in the MEG, the probability of obtain- 
ing more counts than expected in the absorption-line bin 
by chance (assuming the HEG absorption-line equivalent 
width is the mean of a Poisson distribution), is 46.4%. 
Therefore the MEG data do not rule out an absorption 
line. 



probability that the MEG would observe an ab- 
sorption line that is apparently weaker than that 
measured by the HEG. An alternative way of ex- 
pressing the fact that the HEG and MEG data 
are consistent with each other is that the proba- 
bility of measuring between 33-43 counts in the 
absorption-line bin in the MEG is 44.1%. 

We conclude that the significance of the ab- 
sorption line is not low enough to unequivocally 
attribute it to a statistical artifact, but neither 
is the significance high enough to ascribe to an 
astrophysical origin without caution. As we show 
below, the absorption line is significant enough 
to affect modeling of the Fe-K emission line, so 
it cannot be ignored. 

3.2. Gaussian Emission-Line Model 

Here we describe fitting the HEG data (com- 
bined from the +1 and —1 order spectra) with 
a power-law plus Gaussian emission line, with 
and without an absorption line included. If we 
do not include an absorption line, we obtain 
Gaussian Fe-K emission-line parameters entirely 
consistent with those measured by Fang et al. 
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(2002). Namely, we obtained a center energy of 
6.43±g;g| keV, a Gaussian width of 0.10±g;g| keV, 
and an EW of 144+57 e ^ ( au parameters in the 
quasar frame). Note that the EW given by Fang 
et al. (2002) is in the observer's frame so must be 
multiplied by (1 + z) in order to directly compare 
with ours. 

On the other hand, if an inverted Gaussian is 
included to model the absorption line, the emis- 
sion line parameters are different, because there 
is some apparent emission redward of the absorp- 
tion line, that is not modeled when the absorp- 
tion line is not included. This is because a single 
Gaussian model worsens the fit if it is extended 
over the apparent absorption in the data, so it re- 
mains narrow. The power-law plus emission and 
absorption-line model involves a total of eight free 
parameters, and their best-fitting values and sta- 
tistical errors are shown in Table 1. It can be seen 
that, compared with an emission- line only model, 
the peak energy of the emission line shifts down 
by ~ 150 eV, and the intrinsic width is larger, 
because when the an absorption line model is in- 
cluded, there is excess emission on the red side of 
the absorption line that is modeled as part of the 
emission line. Also, since some of the emission 
line is absorbed, a larger intrinsic intensity and 
EW are required. 

Fig. 4 shows the joint two-parameter 68%, 90%, 
and 99% confidence level contours of the ab- 
sorption line EW versus its FWHM. The width 
depends somewhat on how the emission-line is 
modeled (see §3.3 for further discussion). For 
a Gaussian emission line model, the FWHM is 
- 3000 km s" 1 (see Table 1). Fig. 5 (a) shows 
the joint two-parameter 68%, 90%, and 99% con- 
fidence level contours of the Gaussian emission- 
line intensity versus observed center energy, di- 
rectly comparing the two cases, with (black), and 
without (red) an absorption line included in the 
model. The (red) contours, for the model with 
no absorption line, are consistent with those ob- 
tained by Fang et al. (2002), but note that ours 
are plotted against observed energy in Fig. 5 (a) , 
but the same contours are shown against quasar- 
frame energy, in Fig. 9. Fig. 5 (b) shows the EW 
of the emission line versus its FWHM (when the 
absorption line is included in the model). When 
modeled with a Gaussian, the emission line has a 
FWHM of - 20000 km s" 1 (see Table 1). 

We note that Fang et al. (2002) reported the 
detection, albeit marginal, of an additional emis- 
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Figure 4. Joint, two-parameter confidence contours 
(68%, 90%, and 99%) for the equivalent width (EW) 
and FWHM of the high-energy absorption-line feature in 
E 1821+643 when the broad Fe Ka emission line is mod- 
eled in terms of a Gaussian (s olid contours) and in terms 
of relativistic disk line (dotted contours). The absorp- 
tion line is modeled with an inverted Gaussian in both 
cases. In the quasar frame the center energy of the absorp- 
tion line is 6.228^9'?Jg keV (Gaussian emission line) or 

6.220±°;°^ keV (disk emission line). Details of all model 
parameters are given in Table 1 and Table 2. The con- 
tours were obtained from spectral fitting to the HEG data 
(summed over +1 and —1 orders). 



sion line at 6.94_ 07 keV, very likely due to 
Fc xxvi Lya. Using combined HEG and MEG 
data they obtained a AC of 9.2 upon the addi- 
tion of a three-parameter Gaussian to model the 
line, corresponding to a confidence level of 97.3%. 
Fang et al. (2002) measured a line intensity and 
EW of 1.1 ± 0.7 x 10~ 5 photons cm~ 2 s" 1 and 
64 ± 40 eV respectively (we have corrected the 
EW values of Fang et al. for cosmological red- 
shift). Only an upper limit was obtained for the 
intrinsic width, 0.15 keV. Using the same model 
as Fang et al. (i.e. a power law plus two Gaus- 
sians with no absorption line included), we ob- 
tained consistent results for the high-energy iron 
line from the HEG and MEG data. Namely, 
a similar detection significance (98.0%), a cen- 
ter energy of 6.94 ± 0.04 keV, a line intensity 
of O.Stol x 10~ 5 photons cm~ 2 s"\ an EW of 

A0tf 5 eV ' and an intrinsic width of 0.03^;^ keV. 
All statistical errors, ours and those of Fang et al. 
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Figure 5. (a) Contours corresponding to joint two- 
parameter confidence levels of 68%, 90%, and 99% for the 
intensity and center energy (observed frame) of the broad 
Fe Ka line in E 1821+643 when it is modeled with a sim- 
ple Gaussian. Black contours correspond to the case when 
an inverted Gaussian is included to model the absorption 
line in the HEG data (see Fig. 4) and red contours to the 
case when the absorption line is omitted from the model 
and these contours are entirely consistent with Fig. 4 of 
Fang et al. (2002). Caption is abridged: full caption given 
at the end of this article. 



(2002) quoted here for the high-energy iron line 
are 90% confidence for one interesting parameter. 

3.3. Absorption Line Plus Relativistic 
Disk Line Model 

Since the Fe-K line emission appears to 
be asymmetric and peaks below 6.4 keV 
(6.31±g;g| keV, from Table 1), we replaced the 
Gaussian emission line in the model described 
above, with a model of an Fe Ka emission line 
originating in a Keplerian disk rotating around 
a Schwarzschild black hole (e.g. Fabian et al. 
1989). The inner radius of emission was fixed 
at 6r g (r g = GM/c 2 ) and the outer radius was 
a free parameter. The line radial emissivity was 
assumed to be a power law (line intensity pro- 
portional to r~ q ), with the index, q, free. In the 



disk rest-frame the emission-line energy was al- 
lowed to vary between the Fe Ka line energies for 
Fe I and Fe xxvi, allowing a ~ 10 eV margin on 
either side (i.e. 6.39-6.98 keV). The disk inclina- 
tion angle and overall line intensity were also free 
parameters. Thus, there were a total of five free 
parameters for the disk line model. 

The best-fitting parameters of the absorption 
line, disk emission line, and the continuum are 
shown in Table 2, and the best-fitting model is 
shown in Fig. 6 (a). The counts spectrum with 
the best-fitting model folded through the instru- 
ment response and overlaid on the data is shown 
in Fig. 6 (b), whilst Fig. 6 (c) shows the ra- 
tio of the data to the best-fitting model. Note 
that since the C-statistic varied erratically (for 
AC| < 5) for some of the disk line parame- 
ters, Table 2 also shows the 90% confidence, five- 
parameter ranges (i.e. AC = 9.24). The quasar- 
frame absorption-line parameters from the fit are 



6.220 



+0.018 
-0.013 



keV (center energy), 21 



+ 12 



8 eV (Gaus- 
sian width), and 34 ± 13 eV (equivalent width, 
or EW). Fig. 4 shows joint two-parameter confi- 
dence contours for the EW versus FWHM of the 
absorption line (dotted contours), directly com- 
pared with the contours obtained when the emis- 
sion line was modeled with a Gaussian (§3.2). 
The absorption-line parameters are a little dif- 
ferent for the two cases, (notably, the absorption- 
line EW is smaller when the emission is modeled 
with a disk line). 

We compare our best-fitting disk line parame- 
ters to those of Fang et al. (2002). Our inclina- 
tion angle (zero degrees) appears to be smaller, 
but since C varies erratically around this best- 
fitting angle, when one considers the 90%, five- 
parameter upper limit (27 degrees), our result 
is consistent with that of Fang et al. (2002), 

degrees. Our EW of 



who obtained 21.5 jlf 9 4 



209lgy eV is consistent with that obtained by 

Fang et al. (2002), who obtained 168^72 eV > 
which has to be be multiplied by (1 + z) in or- 
der to correct for cosmological redshift. We do, 
however, require a larger line flux in order to 
compensate for the absorption line included in 
the model. We fitted the disk-line model with- 
out the absorption line and obtained a line in- 



tensity of 4.3 



+1.4 
1.4 



x 10 photons cm s 



consistent with the Fang et al. (2002) value of 



3.2 



+ \\ x 10- 5 



photons cm s (again, Fang 



et al. did not correct for cosmological redshift so 
their value has to be increased by 1.297 before 
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Table 2 

Rclativistic Disk Line Fits to Chandra HEG Data for E 1821+643 



Parameter 



Measurement 



C-statistic 

degrees of freedom 

Disk Line Rest Energy (keV) 

Disk Line Emissivity Index, q 

Outer Disk Radius, R out 

Disk Inclination, 6* b s (degrees) 

Disk Line Intensity (10~ 5 photons cm~ 2 s^ 1 ) 

Disk Line EW (eV) 

Absorption Line Center Energy (keV) 
Absorption Line Gaussian Width, a (keV) 
Absorption Line Velocity Width, FWHM (km s _1 ) 
Absorption Line Equivalent Width, EW (eV) 
Power-Law Photon Index, F 
2-10 keV Observed Flux (10~ n ergs cm~ 2 s _1 ) 
2-10 keV Luminosity, Quasar Frame (10 45 ergs s _1 ) 



1014.5 

966 

6 ^7+ 01 
D - ' -o.oi 

(6.51-6.68) 

2 fiq +0.19 

(2.36-3.08) 
> 930 
(> 18) 

0.0183 

(0-27) 
7.0±i"? 

(3.6-10.2) 
209t{& 

(107-305) 

h0.018 
-0.013 
+0.012 



6.220 
0.021_ 

oooc + 1440 
ZJOO -950 
o 4 +13 
°^-13 

1 «4+ - 03 

1 -°^-0.03 

1.2 



3.3 



Simple power-law model plus a relativistic disk emission line, and a Gaussian absorption line, fitted to the Chandra 
HEG data (see §3.3 for details). All parameters (except 2-10 keV flux), refer to the quasar frame (z — 0.297). 
Errors are 90% confidence for one parameter (AG = 2.706); 90% confidence, five-parameter ranges are also given 
for the disk line parameters (i.e. AG = 9.24) in parentheses. Velocities have been rounded to the nearest 5 km s _1 . 
Intrinsic luminosity, in the 2-10 keV band in the quasar frame, calculated using Ho = 70 km s _1 Mpc -1 and 
A = 0.7. 



comparing with ours). 

Although the signal-to-noise of the HEG data 
is limited, the data are quite sensitive to the disk 
inclination angle not being too large and the disk 
outer radius not being too small. If the inclina- 
tion angle is much greater than ~ 30°, and/or the 
outer radius much less than ~ 20r g , then the disk 
line would be too broad to fit the data. We note 
that the inclusion of a Compton-rcflection contin- 
uum to the model described above has a negligi- 
ble effect on the disk line parameters because the 
additional component is not required (|AC| < 1) 
due to poor statistics at the high-energy end of 
the spectrum. For solar abundances, we obtain a 
90% confidence, one-parameter upper limit on the 
solid angle of the reflector of 0.7(27r). For an Fe 
abundance of three times solar, this upper limit is 



0.95(27r). Ionization of the disk and/or inclusion 
of relativistic smearing effects increases this upper 
limit but still has no effect on the line parameters 
because structure in the continuum produced by 
Compton reflection that could potentially affect 
the line parameters, is not required by the data. 

4. Fe-K EMISSION & ABSORPTION IN 
OTHER DATA SETS 

There are two other principal data sets for 
E 1821+643 that we can examine in order to com- 
pare with the Chandra HETGS data, namely data 
from observations with the Chandra low-energy 
grating spectrometer (LETGS) and ASCA. Ex- 
cept for a BBXRT observation, that was very 
short, and in which the Fe-K emission line was 
barely detected (Yaqoob et al. 1993), all observa- 
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Figure 6. (a) Best-fitting power-law plus disk emission- 
line model for the HEG data for E 1821+643, including 
an inverted Gaussian for the absorption line, (b) HEG 
counts spectrum (summed over +1 and —1 orders), with 
the best-fitting model (solid line) in (a) overlaid on the 
data, after folding through the instrument response, (c) 
Ratio of the spectral data to the best-fitting model in (b). 
In both (b) and (c) the spectra are binned at 0.04A. 



tions of E 1821+643 prior to the ASCA observa- 
tion were made with proportional counters, which 
had poor spectral resolution. E 1821+643 has 
not been observed by XMM-Newton (although it 
is present, but too far off-axis, in an observation 
pointed at a nearby brown dwarf). Below we de- 
scribe results from observations with the Chandra 
LETGS and ASCA in turn. 

4.1. Chandra LEG Observation 

E 1821+643 was observed by the Chandra low- 
energy grating spectrometer (LETGS) in 2001, 
January 17-24, resulting in a deep exposure of 
this object. Although the low-energy grating 
(LEG) has a spectral resolution of only 0.05/A 
FWHM (or - 150 eV at 6 keV), it has an effective 
area a factor of ~ 2 larger than that of the HEG at 
the observed energy of the line. The ACIS CCDs 
were used as the detector for the gratings in this 
observation. The low-energy spectrum from this 
observation has been discussed by Mathur et al. 
(2003), but details of the Fe-K complex from this 



data set have not been discussed in the literature. 
We obtained a net exposure time of 470.8 ks, and 
the source showed negligible flux variability dur- 
ing the entire observation. 
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Figure 7. Comparison of the Fe-K complex in 
E 1821+643 between non-contemporaneous Chandra 
HEG, Chandra LEG, and ASCA (SIS) data (see §4). 
Shown are the ratios of data to a simple power-law model. 
Caption is abridged: full caption given at the end of this 
article. 



The ratio of the LEG data (around the Fe-K 
region), to a single power law fitted over the 1-2, 
2.5-9 keV observed energy interval, is shown in 
Fig. 7. An estimate of the line profile itself, with 
the above continuum subtracted off is shown in 
Fig. 8. The latter plot also shows the LEG data 
directly compared with the HEG data, showing 
that the overall line flux in the HEG is somewhat 
larger. It can be seen that there are two apparent 
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emission peaks in the LEG data, not one. One 
peak is centered near 6.4 keV (quasar frame) and 
the other near the energy of the Fe xxvi Lya 
transition, which is at 6.966 keV (rest energy). 
The second emission-line peak in the LEG data 
occurs across several resolution elements and is 
present in both the +1 and —1 orders of the grat- 
ing data. It also has a high statistical significance: 
the difference in the fit statistic between a single 
and double-Gaussian model for the line emission 
is 18.3 (for two additional free parameters in the 
double- Gaussian model). 




4.5 5 5.5 

Observed Energy (keV) 




4.5 5 5.5 

Observed Energy (keV) 

Figure 8. Comparison of the Fe-K complex in 
E 1821+643 between non-contemporaneous Chandra 
HEG, Chandra LEG, and ASCA (SIS) observations (see 
§4). Shown are the absolute fluxes after the best-fitting 
power-law continua (as fitted in Fig. 7) are subtracted. 
Therefore, these plots can be used to assess absolute in- 
tensity changes in the line emission, (a) Chandra LEG 
data (black) compared directly with the Chandra HEG 
data (red), (b) ASCA SIS data (black) compared directly 
with the Chandra HEG data (red). 



LEG data. It can be seen that the energies of 
the emission-line peaks are very well constrained. 
The fact that the first and second emission-line 
peaks in the LEG data are consistent with the 
energies of K lines from Fe I and Fe xxvi respec- 
tively, suggests that they are two separate emis- 
sion lines, rather than a single, broad, disk line 
profile. 

It would appear that the Fe-K emission-line 
profile measured from the LEG is different to 
that measured from the HEG data. However, as 
pointed out in §3.2, Fang et al. (2002) reported 
marginal evidence for a peak near the energy of 
Fe xxvi Lya in the Chandra HEG data, which 
we confirmed (§3.2) and which can be seen in our 
Chandra HEG spectrum (e.g. Fig. 1). In fact, 
formally, the intensity of the high-energy iron 
emission line measured by the HEG is statisti- 
cally consistent with that measured by the LEG 
(compare HEG values in §3.2 with LEG values in 
Table 3). 

No absorption line at the position seen in the 
HEG data is evident in the LEG data. However 
we can obtain an upper limit on the EW of an ab- 
sorption line if it were present, and had the center 
energy and width of the line observed in the HEG 
(see Table 1). We added an inverted Gaussian 
with the HEG parameters to the double-Gaussian 
emission-line model and obtained an upper limit 
on the EW of the absorption line of 23 eV, to 
be compared with the HEG values of 54^^ and 
34^^3 eV, corresponding to Gaussian and disk- 
line models of the Fe-K emission, respectively (all 
values are in the quasar frame). Thus, the LEG 
data do not rule out the presence of the absorp- 
tion line if the disk-line model is more appropri- 
ate for the emission line than the Gaussian model. 
Moreover, strictly speaking, a one-parameter 90% 
upper limit for the LEG absorption-line EW is 
not appropriate since there are a total of eight 
free parameters in the model (plus a continuum 
normalization). The 90%, eight-parameter upper 
limit on the LEG included in the HEG model, 
then the 99% HEG and LEG contours do not 
overlap. Therefore, it is possible that the lower- 
energy peak in the Fe-K line profile varied be- 
tween the HEG and LEG observations. 



Table 3 shows the best-fitting parameters and 
their statistical errors, obtained from fitting a 
power-law plus double-Gaussian model to the 



4.2. ASCA Observation 

E 1821+643 was observed by ASCA in 1993, 
June 19 and the results of this observation have 
been discussed by Yamashita et al. (1997). We 
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reanalyzed the data, in a manner similar to that 
described in Weaver, Gelbord, & Yaqoob (2001). 
We restricted our analysis to the two CCD Solid 
State Imaging Spectrometers (SIS) since they had 
much better spectral resolution than the Gas 
Imaging Spectrometers (GIS) aboard ASCA. At 
6 keV, the spectral resolution of the SIS was 
<~ 150 eV FWHM at the time of the observation. 
We combined data from the two SIS detectors and 
obtained a net exposure time of 83.4 ks. We used 
% 2 statistics for fitting the data since the counts 
per bin are much higher than in the Chandra grat- 
ing spectra, and background has to be subtracted 
for ASCA data. 

The ratio of the SIS data to a single power law 
fitted over the 1-9 keV observed energy interval is 
shown in Fig. 7. This plot also shows the ASCA 
data directly compared with the HEG data. An 
estimate of the line profile itself, with the above 
continuum subtracted off is shown in Fig. 8. The 
latter plot also shows the ASCA data directly 
compared with the HEG data, showing that the 
line appears to peak at a different energy in each 
observation. It can be seen that there is a broad 
emission peak in the ASCA data, centered near 
the He-like Fe triplet transitions (~ 6.6 — 6.7 keV 
in the quasar frame). 

Table 4 shows the results of fitting the ASCA 
data with a power-law plus Gaussian emission- 
line model. The center energy of the line is 
6.62i°;i 2 keV, and the EW is 140t^ eV (both 
in the quasar frame). The fitted to the Chan- 
dra low-energy grating (LEG) data (see Fig. 7, 
Fig. 8, Fig. 9, and §4.1 for details). All param- 
eters (except 2-10 keV flux) refer to the quasar 
frame (z = 0.297). Errors are 90% confidence for 
one parameter (AC = 2.706). The emission-line 
velocity widths (forced to be the same for the two 
lines) were frozen for deriving the error ranges on 
the other parameters, otherwise the fits became 
unstable. The EW of a Gaussian absorption line 
added to the above model is also given, where the 
center energy and width of the absorption line are 
fixed at the respective values obtained from the 
HEG data (see Table 1). Velocities have been 
rounded to the nearest 5 km s _1 . Intrinsic lumi- 
nosity, in the 2-10 keV band in the quasar frame, 
was calculated using H n — 70 km s _1 Mpc -1 and 
A = 0.7. 

Table 4 shows the results of fitting the ASCA 
data with a power-law plus Gaussian emission- 
line model. The center energy of the line is 



6.62i°;^ keV, and the EW is U0±™ eV (both 
in the quasar frame) . The line is unresolved with 
ASCA. Our results are consistent with those of 
Yamashita et al. (1997), who obtained a center 
energy of 6.58 ± 0.05 keV and an EW of 100 ± 
50 eV (both in the quasar frame). Fig. 9 shows 
the 68%, 90%, and 99% contours of line inten- 
sity versus line energy obtained from Gaussian- 
fitting of the emission line in the ASCA data, 
compared with confidence contours obtained from 
the HEG and LEG data. It can be seen that the 
99% confidence contours obtained from modeling 
the principal line peak in the HEG and ASCA 
data overlap (with or without an absorption line 
included in the HEG model). Therefore, at 99% 
confidence, we cannot exclude the case that the 
HEG and ASCA emission-line profiles are con- 
sistent with each other. At 99% confidence, the 
LEG contours for the peak at ~ 6.4 keV also over- 
lap with the ASCA contours. Although the emis- 
sion line at ~ 6.9 keV is not detected in the ASCA 
data, the presence of such a line with the same 
energy and width as the in the LEG data is not 
ruled out by the ASCA data. Adding such an 
additional Gaussian emission line to the ASCA 
data, even the 90%, one-parameter upper limit 
on the intensity, 1.3 x 10~ 5 photons cm~ 2 s^ 1 , 
is consistent with the LEG intensity (0. 8+0.3 x 
10~ 5 photons cm~ 2 s _1 ; see Table 3). 

In order to test the ASCA data for the absorp- 
tion line found in the HEG data, we added an 
inverted Gaussian to our power-law plus single 
Gaussian emission-line model in order to obtain 
an upper limit on the EW of an absorption line 
with energy and width fixed at the values ob- 
tained from the HEG data (see Table 1). There 
were now a total five free parameters, plus a con- 
tinuum normalization. The 90%, one-parameter 
and five-parameter upper limits on the EW are 
82 eV and 138 eV (quasar frame) respectively. 
Thus, the ASCA data do not rule out the pres- 
ence of an absorption line that has the parameters 
that were measured by the Chandra HEG. 

We also modeled the emission line in the ASCA 
data with a relativistic disk-line model, as was 
fitted to the HEG data in §3.3. Since the spec- 
tral resolution and signal-to-noise of the data are 
limited, we fixed r; n = 6r g , r out = 1000r g , and 
q = 2. The line energy in the disk frame, its in- 
tensity, and the disk inclination angle were free 
parameters. We obtained a best-fitting line en- 
ergy of 6.71+0^6 keV (in the disk frame), and an 
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Table 3 

Double Gaussian Line Fits to Chandra LEG Data for E 1821+643 

Parameter Measurement 



C-statistic 
Degrees of freedom 

Ei, First Emission Line Rest Energy (keV) 

11, First Line Intensity (10~ 5 photons cm~ 2 s _1 ) 
EW\, First Emission Line Equivalent Width (eV) 
E 2 , Second Emission Line Rest Energy (keV) 

12, Second Emission Line Intensity (10~ 5 photons cm~ 2 s _1 ) 
EW2, Second Emission Line Equivalent Width (eV) 

a, Gaussian Emission Line Widths (keV) 

FWHM of Emission Lines (km s _1 ) 

Absorption Line Equivalent Width, EW (eV) 

Power-Law Photon Index, T 

2-10 keV Observed Flux (10~ n ergs cm~ 2 s _1 ) 

2-10 keV Luminosity, Quasar Frame (10 45 ergs s _1 ) 



Simple power-law model plus two Gaussian emission lines fitted to the Chandra low-energy grating (LEG) data 
(see Fig. 7, Fig. 8, Fig. 9, and §4.1 for details). All parameters (except 2-10 keV flux) refer to the quasar frame 
(2 = 0.297). Errors are 90% confidence for one parameter (AC = 2.706). The emission- line velocity widths (forced 
to be the same for the two lines) were frozen for deriving the error ranges on the other parameters, otherwise the 
fits became unstable. The EW of a Gaussian absorption line added to the above model is also given, where the 
center energy and width of the absorption line are fixed at the respective values obtained from the HEG data (see 
Table 1). Velocities have been rounded to the nearest 5 km s~ . Intrinsic luminosity, in the 2-10 keV band in 
the quasar frame, was calculated using Ho = 70 km s _1 Mpc -1 and A = 0.7. 
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EW of 16llg7 eV (in the quasar frame), but ob- 
tained only an upper limit on the inclination an- 
gle, of 30° (90%, one-parameter). The inclination 
angle upper limit is consistent with the inclina- 
tion angle inferred from disk-line fits to the HEG 
data (§3.3). However, the ASCA data strongly 
constrain the line energy in the disk frame to be 
near that expected for He-like Fe triplet emission 
(between - 6.6-6.7 keV). The HEG data con- 
strained the line energy to values appropriate for 
Ka transitions from much lower ionization states 
of Fe (i.e. Fe i-Fe xvn or so). This is not sur- 
prising since the HEG data are strongly peaked 
near 6.4 keV, whilst the ASCA data are strongly 
peaked near 6.6 keV (quasar-frame quantities). 

With the above ASCA disk-line model, we ob- 
tained a somewhat smaller upper limit on the EW 
of an absorption line with the parameters mea- 
sured from the HEG data, namely 53 eV (90% 
confidence, one-parameter) . 

An alternative description of the ASCA data is 
that there are two, narrow, unresolved emission 



lines. We fitted a dual-Gaussian model for the 
line emission, with the line widths free but forced 
to be the same value for the two lines. We fixed 
the energy of one line at 6.40 keV (Fe 1 Ka), 
and that of the other line at 6.70 keV (Fe xxv 
resonance, Is 2 — ls2p). We obtained equivalent 
widths of 42I33 eV, and 87lg 4 eV respectively 
(quasar frame). The lower and upper bounds on 
the line intensities from this model are shown in 
Fig. 9, along with intensity versus energy con- 
tours for the Gaussian models fitted to HEG, 
LEG, and ASCA data, as described in previ- 
ous sections. Fig. 10 shows the 68%, 90%, and 
99% joint confidence contours of the intensity of 
one line versus the other from the dual-Gaussian 
model for two unresolved lines in the ASCA data. 
Fig. 10 shows that, in this model, the 6.7 keV 
line dominates the line profile since the 6.4 keV 
is not actually required at 99% confidence and 
higher. We obtained an upper limit of 8,400 
km s^ 1 FWHM on the width of the lines in the 
dual-Gaussian model. 
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Figure 9. Direct comparison of joint, two-parameter, 
68%, and 99% confidence contours of emission-line inten- 
sity, J, versus line energy between non-contemporaneous 
Chandra HEG, Chandra LEG, and ASCA (SIS) data (see 
§4). For ASCA, 90% confidence contours are also shown. 
In each case the data were modeled with a simple power 
law and a Gaussian emission line (two in the case of the 
LEG in order to model the two peaks - see Fig. 7). The 
black HEG contours were obtained when an inverted gaus- 
sian was included to model the absorption feature (the 
same contours are shown in Fig. 5). The red HEG con- 
tours were obtained with no absorption feature included 
(as were the LEG and ASCA contours). The ASCA emis- 
sion line can also be modeled as two discrete, unresolved 
lines at 6.4 keV and 6.7 keV, and the 99% two-parameter 
ranges on the line intensities for this scenario are shown 
as solid vertical bars at those energies (also see Fig. 10). 
Caption is abridged: full caption given at the end of this 
article. 



5. ORIGIN OF THE ABSORPTION 
LINE 

If the absorption line at ~ 6.2 keV detected in 
the HEG data is real, variability is a possible ex- 
planation of non-detection of the absorption line 
in the ASCA and LEG data (but we recall that 
neither of these data sets unequivocally rule out 
the absorption line). If the absorption line is real, 
it could be due to a redshifted resonance transi- 
tion in highly ionized Fe. The high ionization 
state would be consistent with the fact that we 
find no significant soft X-ray absorption features 
in the HEG data (see also Fang et al. 2002), and 
Mathur et al. (2003) came to a similar conclusion 
about the LEG data. If the HEG absorption line 
were due to blueshifted absorption from lighter 
elements, one would expect absorption features 



Figure 10. The ASCA emission line in E 1821+643 can 
be modeled as two discrete, unresolved lines at 6.4 keV 
and 6.7 keV, corresponding to Fe I Ka, Fe XXV Is 2 — ls2p 
resonance transitions respectively (see §4.2 for details). 
Shown are the joint, two-parameter, 68%, 90%, and 99% 
confidence contours of the intensity of the 6.7 keV line 
(^6.7) versus the intensity of the 6.4 keV line (Iba). It 
can seen that the 6.4 keV line is not required at > 90% 
confidence but is allowed to have an intensity high enough 
to be present at the level detected in either the HEG or 
LEG data (see Fig. 9). 



from other ionic species, particularly Fe. If the 
line is due to a Ka transition from Fe xxv or 
Fe XXVI, the observed redshift then corresponds 
to a recession velocity of 20750^^ km s _1 or 
32130iygQ km s _1 respectively (see Yaqoob et al. 
2003b and references therein for the sources of 
atomic data that we use). At the very least, 
v/c ~ 3%, or ~ 9380 km s -1 (corresponding to 
Fe 1). 

At face value, one could interpret the 
absorption- line redshift of ~ 2 — 3 x 10 4 km s _1 
as inflow (e.g. onto the putative central black 
hole). To derive a recession velocity one would 
have to account for gravitational rcdshifts (at the 
unknown radial distances). Evidence for inflow 
is very rare, and in the X-ray band, only one 
other case has been reported (in NGC 3516, Nan- 
dra et al. 1999). In that case inflow was not 
the only interpretation (e.g. see Ruszkowski & 
Fabian 2000) and neither is it for the present case 
of E 1821+643. 

Another interpretation is that the absorption 
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line is due to highly-ionized material that is cross- 
ing the line-of-sight obliquely and may, for ex- 
ample, have have been ejected from the accre- 
tion disk (a similar scenario for line emission was 
considered for Mkn 766 by Pounds et al. 2003a 
and Turner, Kraemer, and Reeves 2004). If the 
ejected blob still has the Kcplcrian motion of the 
disk then from the redshift and upper limit on 
the absorption line width, one can derive a lower 
limit on the radial distance from the black hole 
since, in a given observation time, the blob can- 
not travel too far or else the absorption line would 
be too broad. This gives a constraint on the min- 
imum radial distance of the absorption site from 
the central black hole, as a function of the disk 
inclination angle. This scenario predicts that the 
absorption line energy and EW should be signifi- 
cantly variable. 

Yet another possibility is that the redshift of 
the absorption line is gravitational in origin. The 
absorber could of course be outflowing and yet 
still give an overall redshift, if it is close enough 
to the putative central black hole. From pure 
gravitational redshifts alone, one can place a limit 
on the allowed range of distances, r, from the 
black hole that can be occupied by the absorb- 
ing gas. To first order, from a given line width 
AE, Ar ~ r(r/r g )(AE/E ), where E is the rest 
energy of the line. For r = 10r g , AE ~ 50 eV, 
and Eq ~ 7 keV, we have Ar ~ r g . Obviously, 
the actual allowed range of radii required to keep 
the absorption line narrow depends on the de- 
tails of the geometry of the X-ray source and of 
the kinematics of the absorber (whether inflow or 
outflow). The detailed kinematics of the flow, its 
physical state (e.g. ionization, density, tempera- 
ture), and its proximity to the putative black hole 
will all affect the line width. Nevertheless, Ar is 
rather small. However, we note that the quasar 
model of Elvis (2000) involves a cylindrical out- 
flow (that could intercept the X-ray continuum 
over a narrow range of radii). Reports of out- 
flows in AGN are now common, and a few are 
claimed to originate within tens to hundreds of 
gravitational radii of the central black hole (e.g. 
Pounds et al. 2003b, 2003c). However, it has 
been claimed that many of these outflows may 
not exist if absorption features that are really lo- 
cal to our Galaxy have been misidcntificd as being 
at the redshift of the AGN (McKernan, Yaqoob, 
& Reynolds 2005). 



Table 4 

Gaussian Line Fits to ASCA SIS Data for 
E 1821+643 



Parameter 



Measurement 



Y 2 
A 


1 68 fi 


L»Cgi CCD Ol liet-UOlll 


zuu 


Emission T.iue Rest Enerfv ilceV) 


u - u -0.12 


H.micoinn 1 , i n o mtoyioiT^r 
IJjllllOOlUll J_j1J.1L- lllLUIlolLy 


o 4 +1.9 
^ — 2.1 


(10 — 5 photons cm -2 s — 1 ) 




Emission Line Equivalent Width 


(eV) 


14+ - 16 


cr, Gaussian Emission Line Width 


(keV) 


FWHM of Emission Line (km s _1 ) 


14935±iISlg 


Absorption Line Equivalent Width, 


4±: s 


EW (eV) 


1.80lS;82 


Power-Law Photon Index, F 


2-10 keV Observed Flux 


1.6 


(10~ n ergs cm~ 2 s _1 ) 




2-10 keV Luminosity, 


4.3 


Quasar Frame (10 45 ergs s _1 ) 





Simple power-law model plus Gaussian emission line 
fitted to the ASCA SIS data (see Fig. 7, Fig. 8, Fig. 9, 
and §4.2 for details). All parameters (except 2-10 keV 
flux) refer to the quasar frame (z — 0.297). Errors 
are 90% confidence for one parameter (AC = 2.706). 
The EW of a Gaussian absorption line added to 
the above model is also given, where the center en- 
ergy and width of the absorption line are fixed at 
the respective values obtained from the HEG data 
(see Table 1). Velocities have been rounded to the 
nearest 5 km s _1 . Intrinsic luminosity, in the 2- 
10 keV band in the quasar frame, was calculated 
using H = 70 km s _1 Mpc -1 and A = 0.7. 



6. CONCLUSIONS 

During a Chandra HETGS observation of 
E 1821+643 we detected an absorption line, at 
a significance level of 2 — 3a, at ~ 6.2 keV in 
the quasar frame. Whether or not this absorp- 
tion line is accounted for in modeling the Fc-K 
emission line, directly affects the inferred width 
of the emission line and its peak energy. We have 
also found that the Fe-K emission-line spectra of 
E 1821+643 during observations with the Chan- 
dra HETGS, Chandra LETGS, and ASCA appear 
at first to be different to each other. The first two 
observations were separated by ~ 3 weeks and the 



1G 



T. Yaqoob et al. 



ASCA observation was made nearly eight years 
earlier. The 2-10 keV luminosity of E 1821+643 
during the ASCA and HETGS observations was 
only ~ 60%, and ~ 20% higher, respectively, than 
during the LETGS observation. The Fe-K line 
profile during the ASCA observation was peaked 
at ~ 6.6 keV, during the HETGS observation it 
was strongly peaked at ~ 6.4 keV, and during the 
LETGS observation the line profile was double- 
peaked, at ~ 6.4 keV, and ~ 6.9 keV (all are 
quasar- frame energies). These properties of the 
Fe-K line region are captured in Fig. 7, Fig. 8, and 
Fig. 9. However, due to limited signal-to-noise, 
a non-varying Fe-K complex cannot be ruled at 
99% confidence. It is possible that Ka emission 
from from Fe I-XVII or so, He- like Fe, and Fe XXVI 
Lya is present with similar intensities in all three 
data sets. However, it will be extremely impor- 
tant to test for variability in the Fe-K emission 
complex with future missions. 

Clearly, the structure of the line-emitter must 
be complex (for example, there may be more than 
one emitter), in order to produce this range in 
Fe ionization states simultaneously (e.g. emission 
lines from Fe I- XVII and Fe xxvi are detected with 
a high significance in the LETGS data) . However, 
the location of the line-emitting region or regions 
is not known. Certainly, in the HETGS data the 
Fc-K line at ~ 6.4 keV could be broad, indicating 
an origin close to the putative black hole, such 
as an accretion disk. The spectral resolution of 
the ASCA and LEG data is insufficient to as- 
certain whether the emission lines during these 
observations were narrow or broad, although it 
appears that during the LETGS observation the 
lines were not as broad as the line observed during 
the HETGS observation. 

The putative absorption line detected by the 
high energy grating (HEG) on the Chandra 
HETGS is not unequivocally ruled out by any 
of the data sets. We argue that if the absorp- 
tion line is real, it is most likely due to red- 
shifted resonance absorption in either He-like or 
H-like Fe. The origin of the redshift (equivalent 
to v/c ~ 0.066 — 0.106), could be absorption in 
highly ionized matter that is either inflowing or 
crossing the line-of-sight obliquely. In the latter 
case the material could be due to ejecta from the 
accretion disk. A more controversial interpreta- 
tion is that the absorption line is due to an outflow 
that is so close to the black hole (< 15r s ), that 
the absorption line is gravitationally redshifted. 



If this is true then we would have an impor- 
tant new diagnostic for studying strong gravity, 
complementing studies using the Fe Ka emission 
lines. 

AstroE-2, with a factor of ~ 6 greater effective 
area than the Chandra HEG at the relevant en- 
ergies, and a spectral resolution of <~ 6 eV, will 
help to constrain the origin of the Fe-K emission 
features in E 1821+643, and provide a more sen- 
sitive search for absorption lines of the kind found 
in the HEG data. 
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FULL FIGURE CAPTIONS 
Figure 5 (Full Caption) 

(a) Contours corresponding to joint two- 
parameter confidence levels of 68%, 90%, and 
99% for the intensity and center energy (observed 
frame) of the broad Fe Ka line in E 1821+643 
when it is modeled with a simple Gaussian. Black 
contours correspond to the case when an inverted 
Gaussian is included to model the absorption line 
in the HEG data (see Fig. 4) and red contours to 
the case when the absorption line is omitted from 
the model and these contours are entirely consis- 
tent with Fig. 4 of Fang et al. (2002). The red 
contours are shown again in Fig. 9 of the present 
paper, in that case plotted against quasar-frame 
energy, and these are directly comparable with 
Fig. 4 of Fang et al. (2002). Excluding the ab- 
sorption line results in a slightly higher peak en- 
ergy for the emission line. The vertical dotted line 
corresponds to 6.4 keV in the quasar frame, (b) 
As (a) for the equivalent width (EW) and FWHM 
of the broad Fe Ka line. The dashed line corre- 
sponds to the HEG FWHM spectral resolution at 
the observed peak energy of the broad line, and 
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shows that the Fe Ka line is resolved by the HEG. 
The absorption line is included in the model (see 
Fig. 4). See §3.2 and Table 1 for details. 

Figure 7 (Full Caption) 

Comparison of the Fe-K complex in E 1821+643 
between non-contemporaneous Chandra HEG, 
Chandra LEG, and ASCA (SIS) data (see §4). 
Shown are the ratios of data to a simple power- 
law model, fitted to each data set in the 1-9 keV 
observed energy range (excluding the 2-2.5 keV 
interval for the HEG and LEG data). It can be 
seen that the line profile in each case appears to 
be different, with the LEG profile showing two 
peaks (but see §4.1 and §4.2 for a quantitative 
assessment of variability). The vertical dotted 
lines correspond (from left to right) to the ener- 
gies of the following transitions: Fe I Ka, Fe xxv 
Is 2 — ls2p resonance, and Fe xxvi Lya (that 
correspond to rest-frame energies of 6.400, 6.700, 
and 6.966 keV respectively). It can be seen that 
the HEG profile peaks near Fe I Ka, the LEG 
peaks are near Fe I Ka and Fe xxvi Lya, whilst 
the ASCA peak is near the Fe xxv triplet reso- 
nance line energy. 

Figure 9 (Full Caption) 

Direct comparison of joint, two-parameter, 68%, 
and 99% confidence contours of emission-line 
intensity, I, versus line energy between non- 
contemporaneous Chandra HEG, Chandra LEG, 
and ASCA (SIS) data (see §4). For ASCA, 90% 
confidence contours are also shown. In each case 
the data were modeled with a simple power law 
and a Gaussian emission line (two in the case 
of the LEG in order to model the two peaks - 
see Fig. 7). The black HEG contours were ob- 
tained when an inverted gaussian was included 
to model the absorption feature (the same con- 
tours are shown in Fig. 5). The red HEG con- 
tours were obtained with no absorption feature 
included (as were the LEG and ASCA contours). 
The red HEG contours are consistent with, and 
directly comparable with those shown in Fig. 4 
of Fang et al. (2002). The vertical dotted lines 
correspond (from left to right) to the energies of 
the following transitions: Fe I Ka, and Fe xxv 
Is 2 — ls2p resonance, and Fe xxvi Lya (that 
correspond to rest-frame energies of 6.400, 6.700, 
and 6.966 keV respectively). At 99% confidence 
the evidence for variability in the total line inten- 
sity between LEG and HEG (i.e. for LEG this 



means the sum of the two lines) is marginal, and 
the ASCA 99% contours are large enough to over- 
lap all the other contours. The ASCA emission 
line can also be modeled as two discrete, unre- 
solved lines at 6.4 keV and 6.7 keV, and the 99% 
two-parameter ranges on the line intensities for 
this scenario are shown as solid vertical bars at 
those energies (also see Fig. 10). 
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